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ABSTRACT 

An analytical study was made of the isothermal mass and momentum transfer that 
occurs when a heavy, slow-moving gas is injected coaxially into a duct of light, fast- 
moving gas. The outer stream is composed of two regions,; an outer region moving at 
high velocity, and an intermediate-velocity region adjacent to the inner fluid. As a 
basis for comparison, the total mass flow rate of outer fluid is held constant; thus 
as the buffer region velocity is decreased, the velocity of the remaining outer region 
is increased. The effects of the velocity and thickness of this buffer region on the 
mixing of the two fluids are investigated for both laminar and turbulent flow. 

Outer- t© inner-fluid initial-velocity ratios of 10, 50, and 100 are considered. 
Some representative velocity and concentration fields are shown for laminar and 
turbulent flow. The effect of a momentum buffer region is presented in terms of the 
amount of inner stream fluid that is contained within a specified length of the outer 
duct. The diameter of the duct is taken to be four times that of the jet; duct lengths 
of 2 and 4 jet diameters are considered. 

In a gaseous-fueled nuclear rocket engine, stream mixing is undesirable because it 

tends to dilute the central care of fissionable gas with the surrounding hydrogen 

propellant, and thus leads to an increase in required reactor pressure . The results of 

this study show that both an "optimum" buffer velocity and thickness exist which give 

a maximum amount of inner fluid in the fixed length duct. The optimum velocity is 

greater for turbulent flow than for laminar flow. The optimum buffer thickness is 

shown to be independent of whether the flow is laminar or turbulent. For the system 

studied, the optimum buffer thickness is 1 inner stream radius. The presence of a 
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buffer region increases the amount of inner fluid, in a duct length of 2 jet diameters^ by- 


more than a factor of 2 . 




INTRODUCTION 


The interaction between two coaxially flowing fluids is a subject of considerable 
interest, both academic and practical. A description of mass, momentum, and energy 
transfer between dissimilar coaxial streams is required in order to understand and 
design combustion chambers, jet pumps and ejectors, and afterburners. Recent experi- 
mental and analytical studies have been reported that directly involve such a flow 

pattern. In a supersonic combustor 4 a hydrogen jet is injected into a coflowing 
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stream of oxidizer gas. In a proposed concept for a nuclear rocket engine, a low- 
velocity, fissionable gas is injected into a coaxially flowing, high-velocity stream of 
hydrogen propellant. In both situations, the goal of the studies is to understand, to 
predict, and, ultimately, to control the mixing rate of the two species. 

For the gaseous-fueled-reactor concept, the mixing is undesirable, because it tends 
to dilute the nuclear fuel with hydrogen and thus lead to an increase in required reactor 
pressure. Flow schemes that result in slower stream mixing are therefore of interest. 
Since the mixing occurs primarily because the outer stream velocity greatly exceeds that 
of the inner stream, it is possible, in principle, to delay the interaction of the 
two streams by separating them by means of a buffer region of fluid flowing at some in- 
termediate velocity. This paper describes an analytical study of the characteristics 
of such a flow field. 

In particular, the system to be considered is one in which a jet of low-velocity, 
heavy gas is separated from a surrounding, high-velocity, light gas by an annular sheath 
of intermediate-velocity, light gas. The gases are contained within a duct as illus- 
trative n figure 1. The situation considered is one in which the total mass flow rate 
of the light gas is held constant as it might be in a rocket engine that is required to 
produce a certain amount of thrust. Thus as the velocity of the buffer region is re- 
duced, the velocity of the remaining outer stream must be increased to maintain a con- 
stant mass flow rate. 
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» The analytical method of reference 5 and its extension to turbulent flow 7 have been 

used to investigate the effects of a momentum buffer region on the coaxial mixing of dis- 

► 

similar gases for isothermal flow. The primary parameters are the buffer region thickness 
and velocity. Calculations were carried out for both laminar and turbulent flow. Three 
values of average outer- to inner -stream- velocity ratio were considered. The results are 
compared with the case where no buffer region is present. The results presented indicate 
that a momentum buffer region can be used to alter significantly the rate of mixing in 
the near region of a coaxial jet. 


1,0 


NOMENCLATURE 

b width of mixing region, ft 

c concentration of inner stream species, mole fraction 

I> 1,1 molecular self-diffusion coefficient of species 1, ft^/sec 
d 1,2 binary diffusion coefficient of species 1 and 2, ft 2 /sec 
dimensionless diffusion coefficient 
a constant 

molecular weight, lb/lb-mole 

2r i 0 U 1 0^1 0 

initial Reynolds number of inner stream, 2 1 2 — 

h .,0 

radial coordinate, ft 
dimensionless radius, r/r^ q 
^ q initial Schmidt number of species 1, p^ c/ p l 0^1 1 
axial velocity, ft/sec 
molecular volume ratio 
axial coordinate, ft 
dimensionless axial position, z/r n 
molecular weight parameter, (M-j/Mg) - 1 
eddy diffusivity, ft 2 /sec 
dimensionless eddy viscosity, p e/p 
normalized containment factor 
containment factor 


D 

K 

M 

Re 

r 

r 

Sc 

u 

V 

z 

z 

p 

e 

e + 

^b 
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local viscosity 
dimensionless viscosity 


M- 
P 

p density, lb/ft 3 

\|r stream function (see eqs . 5 and 6) 

Subscripts : 
av average 

b buffer 

max maximum 

w wall 

0 injection point (z = 0) 

1 inner stream species 

2 outer stream species 

i. centerline 

ANALYTICAL METHOD 

The basic analytical method used was reported in reference 5 for laminar flow and is 
extended to turbulent flow in references 6 and 7. It will be briefly reiterated here, 
along with a more detailed discussion of its application to a turbulent three-region 
system. 

Laminar Flow 

The steady-state boundary- layer equations for isothermal axisymmetric flow are: 
Continuity: 

(pvr) + (pur) = 0 (l) 


Momentum : 


du 

v 3F + u 


du 1 d 
5z" - pr Sr: 



(2) 


Diffusion: 


v 



u 


dc p a f rD i p djA 

= r s? V"p 


( 3 ) 


The dimensional equations (l), (2), and (3) are made nondimens ional by normalizing all 
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quantities to a dimensionally similar quantity in the inner stream at the injection 
point. For example, the dimensionless viscosity p is the ratio of local viscosity to 
that of pure inner stream fluid the axial velocity is normalized to the jet 

injection velocity so that u = u/u^q; and - the ax i & l position and radius are normalized 
to the jet radius, z = z/r-|_ q and r = r/r-^Q. The dimensionless diffusion coefficient 
is given by D = (D-^ 2 / D l where D-j^ is the self -diffusion coefficient of the 
inner stream fluid. 


The dimensionless continuity equation is now 


where 3 = (M^/Mg ) - 1. 
the t , z plane, where 
(eq. (4)): 


d_ 

dr 


rv(|3c + l) 


dz 


ru(3c + l) 


= 0 


The system of equations is transformed from the 
the stream function t satisfies the continuity 


(4) 

r, z plane to 
equation 


= -(rv)(3c + l) (5) 

dz 

= (ru)(3c + l) (6) 

dr 

By introducing dimensionless quantities and the stream function along with the 
usual boundary- layer assumption that cty/dz « dty /dr , equations (l), (2), and (3) can be 
written as follows : 

Continuity: 




1 

u(3c + l) 


dt ' 


(7) 


Momentum: 


du _ 2 P + 1 d 
dz R ®i, 0 ^ 


pur (3c + l) 



( 8 ) 


Diffusion: 


dc - 2(3c + 1) _d_ dc 

dz ” Re q ; o^ c l,0 ^ L 


(9) 
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The local viscosity is calculated from a mixing law that can be written in terms of 


dimensionless quantities as 


d = 


Be + 1 


/ \ 1 - c 

(P + l)c + 


( 10 ) 


The dimensionless diffusion coefficient is obtained by writing the Gilliland equation 
( eq . ( 8 ) ) in the form 


r _ D l,2 2 -s/2(P + 2) 


( 11 ) 


where Vg is the ratio of molecular volumes of species 2 to species 1. 

Equations (7) to (ll) are solved numerically by computer program^ to obtain radial 
profiles of dimensionless velocity and concentration at specified axial positions down- 
stream from the injection point. In addition to constants relating to the numerical 
integration scheme, the following inputs are necessary and sufficient to define a case: 

(1) Physical properties: P, dg> Vg, ^ c l,0 

(2) Flow parameter: Re-^Q 

(3) Initial velocity profile: u(0,r) 

(4) Initial concentration profile: c(0,r) 

For all the cases discussed in this paper, flat initial profiles are used, as 
illustrated in figure 1. For concentration, the initial profile is 

0<r<l c = 1 


1 < r < r max 


c = 0 


With a buffer region, the initial velocity profile is 

0 < r < 1 u = 1 


and with no buffer region is 


1 < r < r. 
— b 


r b < r < r max 


0 < r < 1 


1 < r < r 

— max 


U =u b 

u = u. 


2,0 


u = 1 
u = u, 


2,0 
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Turbulent Flow 


The laminar analysis is extended to include turbulent flow by adding turbulent trans- 
port coefficients to their molecular counterparts . Thus the contribution of turbulence 
can be expressed in terms of a turbulent viscosity pe and a turbulent diffusivity e. 

It is also assumed that the eddy diffusivities for mass and momentum transfer are equal. 
The transport properties in equations (2) and (3) are written as 

p t = p + pe (12) 

D i = D + e (13) 

X 

A similar procedure is employed in references 1 to 3, where it is further assumed that 
pe » |i and e » D. 

Utilizing the same mixing law as for equation (10) gives a dimensionless total 
viscosity and diffusion coefficient in terms of a dimensionless eddy viscosity e + 


- Pc + 1 

M t (p + l)c (1 - c) 

1 + e + (1 + e + )p 2 


(14) 


- Pc + 1 

D t (P + l)c + 1 - c 

D + e Sc^ q D + e 0 ( p + 


(15) 


where 


The functional dependence of e + remains to be specified in terms of spatial 
(r,z), flow (Re), and/or property (p,p) characteristics of the system. This has been 
accomplished in current coaxial mixing studies by utilizing Prandtl's postulate that 
in a region of free turbulence the eddy diffusivity is proportional to the width of 
the mixing region and to the velocity difference across it: 

e = Kb ( u max ' 'W < 16 ) 

This postulate and its application to single-component jets and wakes are discussed in 
reference 9. 

The extension of equation (16) to two-component systems can be achieved in 
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different ways. It has been suggested that pe is constant with radiusl>2 and that the 

velocity difference should be replaced with (p 2 u 2 ~ Reference 4 proposes that 

pe is constant with radius and that a combination of mass flux pu and momentum flux 

pu 2 be employed. In reference 10 it is assumed that pe/p. is independent of radius, 

but dependent on axial position to some exponent. Reference 7 concludes that pe/p 

can be taken as independent of radial and axial position and that the viscosity 

l/S 

difference should be expressed as ( | u s/ u l " "H ' * 

Although these various approaches have yielded unresolved differences, they all 
share some significant similarities. All the suggested formulations (l) are in 
general agreement with experimental data, (2) are applicable for initial velocity 
ratios above and below 1, (3) indicate that the radial dependence of pe is negli- 
gible, and (4) conclude that density ratio Pg/pp must be included as a parameter 
when applying equation ( 16 ) to two-component systems. These similarities suggest 
that each approximation will correlate data within certain ranges of flow and property 
parameters and that the difference between laminar and turbulent flow is much greater 
than the differences between the various algebraic descriptions of turbulent viscosity. 
The numerical solution of reference 5 and the turbulent viscosity expression of refer- 
ence 7 have been used for the calculations reported herein because no linearization 
assumptions are required and because general agreement with turbulent coaxial mixing 
data has been demonstrated over the widest range of initial velocity ratios 
(0.83 < u 2 <0.97 ; 1.25 < u 2 < 49) . 

A recent study has shown that these various correlations predict about the sm 
turbulent viscosity at the jet origin when they are compared on a consistent basis. 

It is also shown that, for nearly equal stream velocities, turbulent mixing can be 
reduced by the presence of honeycomb inserts to remove "per turbulence" from the two 
streams just prior to the jet origin. The turbulent viscosity correlation of 
reference 7 used here included both the perturbulence and that due to the shear between 
the two free streams . 

From reference 7, the ratio of turbulent to laminar viscosity can be written as 
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(17) 


e+s f= °- i2 (^XD 8ei ’° - 25 ° )(u 2 - 1,1/2 

* For a two-region coaxial-flow situation, ug is simply the initial velocity ratio 
of the outer to the inner fluid, and the use of equation (17) is straightforward. When 
a momentum buffer is present, there are two velocities in the outer fluid, and it is not 
so clear how to obtain one velocity to represent all of the outer fluid. Certainly this 
velocity should be some kind of average of the buffer region velocity and the outer 
region velocity, but there are a number of possibilities: a simple numerical average, 

a mass -flow-rate-weighted average, and a momentum-weighted average. The latter two 
procedures seem more realistic than a numerical average. For the most extreme case 
considered, the buffer region contained only 10 percent of total mass flow of the 
outer fluid, and 4 percent of the total momentum. Thus either of the two weighted 
averaging methods produced a velocity that was essentially the same as that of the 
high-velocity outer region. The choice, then, is rather academic; the u 2 in equa- 
tion (17) was evaluated with a momentum -weighted average velocity of the outer fluid. 

SCHEDULE OF CALCULATIONS 

The primary variables to be considered are the thickness of the buffer region, 
the velocity in the buffer region, and the initial velocity ratio of light gas to 
heavy gas. All calculations were made for both laminar and turbulent flow. 

The first of each series of cases was for no buffer region present. For example, 
an initial velocity ratio of 10 was selected for the first case. For the next case, 
a buffer region thickness of 0.5 initial inner stream radii (r^ = 1.5) was assigned, 
and the velocity ratio in this region was decreased until it was some specified 
fraction of that of the remaining outer fluid. A dimensionless channel radius r w of 
4 was assigned for all calculations, and as the buffer velocity was decreased, the 
velocity of the remaining outer fluid was increased to maintain a constant total mass 
flow rate of light gas. Thus if the initial velocity ratio with no buffer region is 
u 2 , 0 , then 

"2,o(*v - ^ = ^ + "2,0<*v ' *b) (1Q ) 
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For the illustrative case selected, u 2 ^ 0 = 1° and r b ~ 1*5; for a buffer velocity that 
is 0.1 of the remaining outer fluid velocity. 



and u-^ = 1.08. 

Average outer- to inner -velocity ratios of 10, 50, and 100 were considered. For 
each of these ratios, buffer- to outer-velocity ratios of 0,1, 0,2, 0,3, and 0,4 were 
investigated. For all combinations of these variables, buffer thicknesses of 0.5 and 1 
were calculated, as well as the situation of no buffer region. These combinations were 
calculated for both laminar and turbulent flow: 

Average outer- to inner-velocity ratio, Ug^Q, 10, 50, 100 
Buffer- to outer -velocity ratio, u^/ug^o* 0.1, 0.2, 0.3, 0.4 
Buffer region radius, r^, 1, 15., 2.0 
Flow regime, laminar, turbulent 

For all calculations, the following constant values were used: 0 = 100, jig = 0.2, 

Vg = 0.2, Sc-^ q = 1.0. As pointed out in references 5 and 10, the analytical results 
are relatively insensitive to the physical-property ratios, so these choices are not 
critical. Some of the calculations performed for this study were repeated for molecular 
Schmidt numbers of 0.5 and 2, and the concentration and velocity profiles were virtually 
unaffected. Insensitivity of a similar analysis to Schmidt number was reported in 
reference 4. The laminar calculations are based on an inner stream Reynolds number of 
200 and the turbulent cases for a Reynolds number of 20 000. 

RESULTS AND DISCUSSION 

The results of the calculations of three-region, two-component mixing will be 
presented in two main categories. First, detailed characteristics of concentration 
and velocity fields will be discussed for both laminar and turbulent flow. Next, a 
more generalized presentation of the effects of a buffer region will be made in terms 
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of the total mass of inner stream fluid present in a channel of specified length. 

► 

Flow Field Characteristics 

The general effect of a faster moving outer stream is to accelerate and dilute 
the inner stream fluid. This results from mass diffusion of the inner fluid radially 
outward into regions of higher velocity, and also from acceleration of the inner fluid 
by momentum transfer radially inward. Figure 2(a) shows the concentration field for 
turbulent coaxial flow with an initial velocity ratio of 50. If the initial velocity 
ratio were lower than 50, the mixing of the two streams would occur more slowly. A 
buffer region of lower velocity immediately adjacent to the inner stream presents such 
a situation, at least initially. 

Figure 2(b) shows the effect of reducing the velocity in a buffer region, 1 jet 
radius in thickness, until it is 0.1 that of the remaining outer stream. In order to 
maintain the same mass flow of outer fluid, or the same average velocity ratio of 50, 
the outer-fluid-velocity ratio between the buffer region and the channel r w =4 is 
increased to 61 and the buffer velocity ratio is 6.1. Figure 2(b) illustrates that the 
general effect of a momentum buffer is to reduce the mixing rate in the near region 
of the jet. Farther downstream, the mixing zone produced by the shear between the 
high velocity outer hydrogen flow and slow moving buffer region penetrates to the 
inner zone and results in increased acceleration of the inner fluid. 

Figure 2(c) shows the concentration field that results if the buffer-to outer- 
velocity ratio is increased from 0.1 to 0.3. Although the buffer velocity is higher, 
the mixing rate is less than for a buffer-velocity ratio of 0.1. This indicates that 
there is an optimum buffer velocity and that buffer velocities below this value do 
not afford sufficient momentum in the buffer to offset the increased momentum in the 
outer stream required to maintain a constant mass flow. A comparison of the concen- 
tration lines c = 0.6 in figures 2(a) to (c) shows this effect. Figure 3 is a plot 
of centerline concentration as a function of axial position for the three cases first 
discussed. The centerline concentration is higher for a buffer -velocity ratio of 0.1 
for only l/2 jet radius downstream; beyond 2 radii downstream it falls below that for 
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no buffer region. Between axial distances of 1 and 4.5 jet radii, a buffer -velocity 
ratio of 0.3 affords more protection of the inner stream. Figure 3 also serves to illus- 
trate that, although centerline concentrations exhibit an effect of a buffer, it is 
rather minimal. The effect of a buffer region on the velocity field is quite pronounced. 
This is shown in figure 4, again for turbulent flow and an average outer- to inner- 
velocity ratio of 50. Figure 4(a) is for no buffer region, and figure 4(b) is for a 
buffer thickness of 1 jet radius at a buffer- to outer-velocity ratio of 0.3. The 
reduction in the rate at which the inner fluid is accelerated is clearly shown by a 
comparison of the cons Lant velocity lines for u = 5 (fig. 4). The constant velocity 
lines for u = 49 show that, farther downstream, the higher velocity outer fluid is 
beginning to overcome the initial buffer effect. 

Figure 4 shows the entire velocity field regardless of species and does not 
threfore clearly reflect what has happened to the inner stream fluid. Figure 5 shows 
axial variation of the radially averaged velocity of the inner fluid, still for turbu- 
leng flow, a buffer thickness of 1 jet radius, and an average initial velocity ratio 
of 50. The effect of a buffer region is shown more clearly in figure 5 than the center- 
line concentrations showed in figure 3. A buffer- to outer -velocity ratio of 0.3 has 
reduced the rate of acceleration of the inner stream fluid. At a distance of 8 jet 
radii downstream, the presence of a buffer with a velocity ratio of 0.1 has resulted 
in a higher velocity inner fluid than the case with no buffer region. 

The influence of a buffer region on the average inner fluid velocity can be 
interpreted as an effect on the mean residence time of the inner fluid in a given length 
of channel. Thus the effect of a buffer region can be expressed in terms of its effect 
on the total mass of inner fluid that is contained in a channel of length z, compared 
with the maximum amount that could be present. This maximum would occur for slip flow 
and no radial diffusion. The ratio of the total amount of inner fluid present to that 
which would exist with no interaction between the two streams is then a sort of "contain 
ment efficiency", or containment factor. It represents the degree to which the inner 
stream fluid has been swept from the system. It can vary from a maximum of 1.0 down to 
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l/ug^ 0-> which would represent instantaneous acceleration of the inner stream fluid. 

Figure 6 shows a plot of this containment factor T] c for various channel lengths 
z , turbulent flow, a buffer thickness of 1 jet radius, and u/g q = 50. Figure 6 
shows the significant effect of a buffer region on the amount of inner stream fluid 
contained in a channel of length z . For a system that is 4 jet radii in length, a 
buffer region with a buffer- to outer -velocity ratio of 0.3 increases the amount of 
inner fluid present from 18 percent without a buffer to 31 percent of the theoretical 
maximum. A buffer -velocity ratio of 0.1 yields a containment factor of 0.26. 

Figures 7 and 8 show turbulent concentration fields for initial velocity ratios 
of 10 and 100, respectively. Comparison of these figures with figure 2 (for Ug^o = 50) 
indicates that the effect of a momentum buffer is more pronounced at higher initial 
velocity ratios, quantatively shown in succeeding figures. 

For an initial velocity ratio of 50, the ratio of turbulent to laminar viscosity is 
830 with no buffer region; this value was used in the calculation illustrated in fig- 
ure 2(a). With a buffer region present, a momentum-averaged velocity of the outer fluid 
yields a value of 917 for pe/p for the case shown in figure 2(b). Figures 9(a) and (b) 
show the corresponding concentration fields for laminar flow (pe/p = 0) for these two 
cases. A comparison of the concentration lines for c = 0.2 in figures 9(a) and (b) 
shows that a momentum buffer has a considerable effect for laminar flow, even though the 
general nature of the flow fields for laminar and turbulent mixing is markedly different. 

Effect of Buffer Parameters on Containment Factors 

Containment factors for laminar and turbulent flow are shown in figure 10 for basic 
coaxial mixing, with no buffer region present, and for initial velocity ratios of 10, 

50, and 100. As before, the containment factor is the ratio of the total mass of jet 
fluid contained in a duct of length z to the amount that would be present with no 
interaction between the two streams . This containment factor serves as a measure of the 
degree of stream mixing. Higher initial velocity ratios, longer ducts, or the presence 
of turbulence all lead to increased mixing and lower containment factors. 

The presence of a momentum buffer region reduces stream mi xing and increases the 
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containment factor. The effects of buffer-region thicknesses and velocities on the con- 

** 

tainment factor are shown in figure 11 for laminar flow and in figure 12 for turbulent 
flow. Curves are shown for buffer-region thicknesses of 0.5 and 1.0 jet radii and for 
initial average velocity ratios of 10, 50, and 100. 

The curves exhibit the same general effect of buffer velocity on containment factor 
for laminar and turbulent flow. As the buffer-region velocity is decreased, the contain- 
ment factor increases to a maximum value and then begins to decrease. A comparison of 
two corresponding curves for z values of 4 and 8 shows that the buffer effect is re- 
duced for longer ducts, but that the optimum buffer- to outer-velocity ratio is essen- 
tially the same. The remaining results will be presented for a duct that is 4 jet 
radii in length (z = 4). 

It is difficult to assess the relative effects of a buffer region from the various 
curves of figures 11 and 12 because the absolute value of the containment factor is 
strongly affected by the initial velocity ratio and by the nature of the flow. It is 
therefore useful to normalize each curve to its value with no buffer region 

(*b/u 2 ,0 = 1-0). 

The variation of such a normalized containment factor ti, with buffer- to outer- 

'b 

velocity ratio is shown in figure 13. The parameter rj-^ is the factor by which the 
amount of jet fluid in the duct has been increased by a buffer region. For an average 
velocity ratio of 100, a buffer that is 1 jet radius thick (r^ = 2) and for turbulent 
flow, the amount of jet fluid in the duct is increased by a factor of 2.25 by decreasing 
the buffer-velocity ratio u^/u^ ^ from 1 (no buffer) down to 0.2. The maximum buffer 
effect is obtained for a buffer-velocity ratio of approximately 0.2, though there is 
some variation. Figure 14 shows the optimum buffer -velocity ratio as a function of the 
average initial velocity ratio Ug q. For higher initial velocity ratios, the optimum 
buffer-velocity ratio is lower. In general, the optimum buffer velocities range from 
0.2 to 0.4 of the outer velocity for turbulent flow and from 0.1 to 0.3 for laminar 
flow. Thicker buffer regions require lower buffer velocities. 

Figure 15 shows the effect of the thickness of the buffer region on containment 


14 



factors that have been maximized with respect to the buffer-velocity ratio. As with 

buffer-velocity ratio, an optimum exists. Higher initial average velocity ratios require 
* 

somewhat thickner buffer regions, but the effect is relatively weak. For the range of 
conditions studied, the optimum momentum buffer thickness is approximately 1 jet radius. 
-It is interesting to note that the optimum buffer thickness is virtually independent 
of whether the flow is laminar or turbulent. 

CONCLUSIONS 

The results of this study indicate some characteristics of two-component, three- 
region coaxial flow. All the calculations were carried out for a channel- to jet- 
radius ratio of 4, and for the following physical properties : A ratio of inner fluid 

to outer fluid molecular weight a ratio of 100, an inner fluid Schmidt number of 1, 
a ratio of outer-fluid to inner-fluid viscosity of 0.2, and a ratio of outer-fluid 
to inner-fluid molecular volume of 0.2. For laminar-flow calculations, the jet Reynolds 
number was taken as 200; for turbulent flow, the jet Reynolds number was taken to be 
20 000. Channel lengths of 4 and 8 jet radii were considered, and buffer thicknesses 
of 0, 0.5, and 1 jet radii were investigated. Average ratios of outer-fluid to inner- 
fluid initial velocity of 10, 50, and 100 were studied; for each of these cases, buffer- 
to outer-initial velocity ratios of 0.1, 0.2, 0.3, and 0.4 were considered. For these 
ranges of conditions, the following results were obtained: 

1. An optimum buffer-region velocity exists that minimizes the interaction between 
the two fluids , 

2 . Similarly, there is an optimum buffer-region thickness . 

3. The optimum buffer velocities range from 0.2 to 0.4 of the outer velocity for 
turbulent flow and from 0.1 to 0.3 for laminar flow. 

4. The effect of a momentum buffer region is greater for higher initial velocity 
ratios of outer stream to inner stream fluid. 

5. For a duct radius that is four times that of the jet, the optimum buffer thick- 
ness is 1 jet radius for both laminar and turbulent flow. 

6. The presence of an intermediate -velocity momentum-buffer region of outer fluid 
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significantly affects mass and momentum transfer and can increase the amount of inner 
fluid contained in a duct 4 jet radii in length hy more than a factor of 2. 
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u 2 , 0 = 50 * 



Y/////////////////, 

(a) No buffer region. 


u 2 0 = 61 * 



(b) Buffer thickness, 1; q = 0- L 
Figure 2. - Effect of buffer region on concentration field. Turbulent flow; 
520 * 5 °. 
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(c) Buffer thickness, 1 ; Ujj/u 2 q = 0. 3. 

Figure 2. - Concluded. 



Figure 3. - Centerline concentration. Turbulent flow; 
u 2 0 - 50; buffer thickness, 1. 








Axial position {initial inner radii), z 


Figure 5. - Effect of buffer region on inner fluid average 
velocity. Turbulent flow; u 2 g - 50; buffer thickness, 1. 



Figure 6. * Effect of buffer region on total mass of inner 
fluid in system of length z. Turbulent flow; u 2 g * 
buffer thickness, 1. 
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(a) Axial position, z/r^ q ■ 4 ; aver- 
age initial velocity ratio, U 2 q * 10. 



(c) Axial position, z/r 1 q = 4; aver- 
age initial velocity ratio, u 2 q * 50. 

4r- 


Buffer to jet- 
radius ratio, 

? b 



(b) Axial position, z/rj 0 * 8; aver- 
age initial velocity ratio, S 2 o * 10 * 



(d) Axial position, z/r^ q = 8; aver- 
age initial velocity ratio, u 2 0 * 50. 



0 .4 .8 1.2 0 _.4_ .8 1.2 

Buffer- to outer-velocity ratio, u^u 2 q 


(e) Axial position, z/rj 0 * 4; aver- if) Axial position, z/r 1 0 = 8; aver- 
age initial velocity ratio, u 2 q b 100. age initial velocity ratio, u 2 0 - 100. 

Figure 12. - Buffer containment factors. Turbulent flow. 







Buffer- to outer-velocity ratio, u^/ u 2 q 
< b) Buffer- to jet-radius ratio, r b - 2.0. 


Figure 13. - Normalized buffer containment factors for 
z ,r i 0 -4. 
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Figure 14. - Optimum-buffer to primary-velocity ratios for z/rj 0 - 4. 



Figure 15. - Optimum buffer thicknesses for 

z,r 10*4. 
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